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The development of Lie methods to optics is very recent. for a 
review of this approach see [1]. Through studying linear and 


nonlinear transformations (aberrations) of optical phase space 
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which mode! optical systems, Lie theory can provide an effective 
calculation method in optics [2,3,4]. Light rays are described in 
an optical phase space as points (Pie Pre gdp) evolving along the 
optical axis z of the system ( which takes the role of time in the 
hamiltonian formulation of classical mechanics). At every 
z=constant plane the configuration space has coordinates 
aatar aD Fermat's principle leads to the optical Lagrangian and 
allows to calculate the conjugate momenta Pairs which can he 
interpreted as a two-dimensional vector p in z=const. plane, 
along the projection of the ray on the plane. Its magnitude p is 
determined by the refraction index n at (q,z) and the angle @ 
between the ray and z axis: pensing. A general linear group which 
action conserves the optical phase space structure is a symplectic 
group Sp(4,R). For axially symmetric systems it is enough to use a 
smaller symplectic group, namely Sp(2,R). However. for a general 
case of nonmaxial optical media one has to work with full Sp(4,R) 
group which has rather complicated structure. Instead of Sp(2,R) 
as a subgroup of Sp(4,R) it is more convenient, in cas case, to 
use a chain A a maximal compact subgroup, a very well 
known unitary group U(2) c Sp(4.R). Embeding of this group into 
symplectic one is given by the construction [5]: 


if A+iB € U(2), where A and B are two-dimengional real matrices. 
AB 3 ; 
then -p al © U(2) c Sp(4,R). (1) 


The generators of the group U(2) can be calculated directly from 


(1) in the form of 4x4 matrices: 
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where GA k=1,2,3, are the usual Pauli matrices and Po denotes the 
2x2 identity matrix. A subgroup U(1) < U(2) is generated by the 
matrix Xa Three matrices Xp’ Xe and Xa generate the subgroup 
SU(2) . Two diagonal generators 


D,diag(1,0,-1,0) and D.=*diag(0,1.0,-1) (3) 


2 
„can be easily found from Iwasawa decomposition of Sp(4,R) [6]. 
Last four generators can be calculated from the following 
comutators: 

X,-31X,-Dy]- X--51X, Dol 

X57(X,-Do} > Xa” (Xa D1) (4) 
This way we have constructed required generators for Sp(4,R)> U(2) 
group chain. 
First approximation to optical ray dynamics is a linear optics. 
We consider here an optical medium described by the quadratic 
refraction .index 


nen + DA gaa (5) 
o u i j 5 


where Na’ we assume to be negative numbers, to obtain an 
“atractive" optical system. For the oposite case one can calculate 
the evolution operator analogously. The optical hamiltonian for 
linear optics can be thus written. in the form 
UL * 
H 5a P n>) naa . Ke? 
o tj 

Lie operator corresponding to a function on the optical phase 
space we denote by ~ , e.g. for H Lie operator will be denoted as 


H. In our case H is a simple differential operator acting on 


functions defined on the optical phase space: 
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rio aaa 3p, t aael Ha 
Action of H onto column vector (P) -P3949393 allows to find its 
matrix representation, which in turn can be expressed in terms’ of 
the generators (2-4): 
H = 4(n,,+noo-n*)X,+ $(n,.-n55)X,+ $in, +n) X 
Bee gS AI 22 oita ee ae Dome et 
Sean, tn) X,+ $n X fnat) Xg (8) 
The hamiltonian H and its representations (7) and (8) are 
independent of z ("optical time"). This property implies that the 
evolution operator Ga (Zz) = exp(-z H ) ( in the matrix 
represetation). The operator Gyu (2) can be calculated directly 


making use of Perron's formula for n-th power of matrix E : 


eS i aede = ji F(A) ] By 
= {a@,-I)! ail AEF 
os at da ROOS AmA, 
i i 
where A (i=1,2,...,t) denotes different eigenvalues of the. matrix 


EN, a their multiplicity and F(A) is an algebraic complement of 
the matrix AIE. The evolution operator can be written in a form 


of 4x4 z-dependent matrix - 


"E G, (2) Wie é (10) 
where (Wy W2. Wg- W4) =(P] -P393 +49) - For the refraction index (5) 


with negative Nyy and Noo (only this case is interesting for light 
Propagation in optical, fibers ) one can get, in general, four 
eigenvalues `j which are puerly imaginary numbers: 

ayr (DE Tig ( (xt (1-265, 25a) v1 /n 2, an 
where x=-(n,,+nN55)>0 and Yelm a a ode. « In this 
case the matrix elements of the evolution operator can be 


expressed in terms of trigonometrical functions 
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G |, (2) = Ay, cos (94 Z) +B, cos(w5z) +O, Sin (wz) +D Sintez). (12) 


where the spacial frequences are given by a" L(x (-1) *yy /n 11? 


and all non-zero coefficients we list below: 
A, 47Ag3770-S (x-y+2n55)/y By 17B3370-5(x+yY+2n33)/Y 


Agg™Agg™—0-S(x-yt2n,,)/¥ B357B =0.5(x+y+2n,4)/¥ 


44 
Bio gi ag hase eo ates aes 4am 2 200 
Czo"Cg,7-0.5 (ny 3+3) /tY in, (x-y) 1**) 
Dyo7D 4170-5 (ny +44) /{yin, (x+y) 1} 
Cy g2Co570-5(n,o#n5,) [n (x-y) 1°77 
Dy 470237-0 -5 (ny otno,) In, (x+y) 1° /y 
CaO. 5n 7 lAn] naa (n4 2+3) ten, OY) Vly (xsvie lk 


} 


4/2 


Dysr-0-Sn JE 
3 a 


[4n, n327 (ny o+n,) +20} (x+y) 1/Ly (x+y) 
C3470. 5ni [4n, poo (0124031) +2033 (xy) 1/ Ly (x-y) 7] 
Da47-0.505 7 lAn; yng 5-(ny atng) 2ng (x+y) 1/ Ly (xty) 7] 
C3170.5(21n33+x-Y)/ {ylin (x-y) 17) 
D4, 7-0,5(2n55+x+y) / tyin (x+y) 1°77) 
Cyo70.5(2n, ,+x-¥) /(y tn, (x-y) 1) 
Dao -0-5(2n,,+x+v) / (yin, (x+y) 1°77} 
From the expression (12) one can see that motion of light rays is 
much more complicated in this case than in presence of axial 
symmetry. However, even now the motion without aberrations is 
harmonic in the optical phase space. To include aberrations one 
can follows the works [2,3] and [4]. Using the basis constructed 
in [4] and the algebra a? non-linear effects in light ray motion 
can be approximately taken into account. The formula (12) gives a 
mean part of light ray motion and can be useful in preparing of 


more exact numerical program. Inclusion of third order aberrations 
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allows to obtain more. details about light rays propagation in 
quadratic optical media, especially in optical fibers. Solution of 
this problem requires a construction of appropriate 
representations of the optical group Sp(4,R) and will be presented 


in subsequent paper. 
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